Superfluid Turbulence: Nonthermal Fixed Point in an Ultracold Bose Gas 
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Nonthermal fixed points of far-from-equilibrium dynamics of a dilute degenerate Bose gas are 
analysed in two and three spatial dimensions. For such systems, universal power-law distributions, 
previously found within a nonperturbative quantum-field theoretic approach, are shown to be re- 
lated to vortical dynamics and superfluid turbulence. The results imply an interpretation of the 
momentum scaling at the nonthermal fixed points in terms of independent vortex excitations of 
the superfluid. Long-wavelength acoustic excitations on the top of these are found to follow a non- 
thermal power law. The results shed light on fundamental aspects of superfluid turbulence and 
have strong potential implications for related phenomena studied, e.g., in early- universe inflation or 
quark-gluon plasma dynamics. 
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From the formation of Bose-Einstein condensates in ul- 
tracold gases to quark-gluon plasmas produced in heavy- 
ion collisions, over a range of twentyfour orders of mag- 
nitude in temperature, nonequilibrium dynamics governs 
many interesting phenomena. Turbulence is an outstand- 
ing and intricate example, which can be described as an 
oriented stationary flow of energy or particles between 
different scales. This idea of an energy cascade goes 
back to the work of Richardson in the context of at- 
mospheric science Kolmogorov, in his 1941 math- 
ematical discussion of turbulence in an incompressible 
fluid, added the concept of universality and scaling 0. 
More recently, turbulence has been studied in the context 
of the inflationary early universe as well as of strongly 
correlated matter produced in heavy-ion collisions 13|- 
01 -Dynamical critical phenomena, i.e., fixed points of 
the evolution away from thermal equilibrium have been 
proposed. They potentially affect the equilibration pro- 
cess by forcing the evolution to critically slow down be- 
fore final thermalization. New scaling laws were found 
by analysing non-perturbative Kadanoff-Baym dynamic 
equations 0,0 ■ Analogous predictions for a nonrelativis- 
tic Bose gas were given in [6[, proposing strong matter- 
wave turbulence in the regime of long-range excitations. 

Superfluid turbulence, also referred to as quantum tur- 
bulence (QT) has been the subject of extensive studies 
in the context of helium Q . In contrast to eddies in clas- 
sical fluids vorticity in a superfluid is quantized [lol [Tlj , 
and the creation and annihilationprocesses of quantized 
vortices are distinctly different The observation of 
a Kolmogorov 5/3-law 0] in experiments with superfluid 
helium, cf. (l^ for a review, received much attention [Tsl - 
[Toj . In particular, the role of the normal-fluid as com- 
pared to the superfluid component in the turbulent flow 
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is under debate Q . In the context of the kinetics of con- 
densation and the development of long-range order in a 
dilute Bose gas, the role of turbulence in the su per fluid 
and its acoustic excitations was discussed in Refs. |2Cll . [2l| . 
see also [22| for a recent review. A possible observation of 
QT in ultracold atom ic g ases presently poses an exciting 
task for experiments |23l [24|. 

A satisfactory ab-initio mathematical description of 
both quantum and classical turbulence is inherently dif- 
ficult due to the generically strong correlations building 
up within the system. Analytical results are known, how- 
ever, in regimes where kinetic theory applies: In a dilute, 
degenerate Bose gas the normal-fluid component can vary 
at the expense or gain of the superfluid part. As a conse- 
quence, the gas is compressible and so-called weak wave- 
turbulence phenomena can occur. For these, scaling laws 
can be derived by analysing kinetic equations j25j . 

Recent developments presented in Refs. 13 allow 
to set up a unifying description of scaling, both in the ul- 
traviolet quantum-Boltzmann kinetic regime and in the 
infrared (IR) regime of long wavelengths where perturba- 
tive approximations break down. In the IR regime, new 
scaling laws were found Q by analysing non-perturbative 
Kadanoff-Baym dynamic equations with respect to non- 
thermal stationary solutions. 

Here, we study QT in two and three dimensional dilute 
Bose gases by means of simulations in the classical-wave 
limit of the underlying quantum field theory. We 
compare the scaling exponents of the single-particle mo- 
mentum distribution and their analytical predictions . 
While we find excellent agreement, our results provide 
an interpretation of the dynamical fixed points proposed 
in Refs. [1, Q for the case of an ultracold Bose gas: The 
appearance of nonperturbative infrared scaling reflects 
the presence of independent vortices. Moreover, the 
power spectra of underlying compressible excitations 
suggest an understanding in terms of acoustic turbulence 
[251 on the top of the vorticity-bearing quasicondensate. 
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Scaling laws: Analytical Predictions. Based on the con- 
cept of universality, the turbulent single-particle spec- 
trum 71 (k) is assumed to scale as 



7i(sk) 



n(k), 



(1) 



where s is some positive real number. This scaling law 
is compatible with a more general scaling behaviour as- 
sumed for the full two-point correlation function of the 
bosonic quantum field as discussed in Ref. @. 

To determine the positive exponent C in the infrared, 
where n(k) ~ |k|^^ is large, an approach beyond ki- 
netic weak-wave-turbulence theory has been used. This 
is based on analysing fixed points of Kadanoff-Baym dy- 
namic equations derived from the two-particle irreducible 
(2PI) effective action or ^-functional [2^. While in the 
regime of large wave numbers the approach goes over 
into the kinetic description of weak wave turbulence, 
an effectively renormalised many-body T-matrix modi- 
fies the scaling law in the IR. This is a consequence of 
the resummation of a particular class of Feynman di- 
agrams contributing to the 2PI effective action [1, Q , a 
formulation substantially beyond the low-order perturba- 
tive expansion used for the kinetic quantum Boltzmann 
equation. Physically, the renormalised T-matrix implies 
a reduction of the effective interaction strength in the 
IR regime of strongly occupied modes As a conse- 
quence, single-particle occupation numbers rise, towards 
smaller wave numbers, in a steeper way than in the weak- 
turbulence regime. The IR scaling exponent for constant 
radial quasiparticle flow in d dimensions was predicted in 
Ref. [I as 



(2) 



Semiclassical field simulations. A dilute superfluid 
Bose gas can be described, in the classical-wave limit, 
by the Gross-Pitaevskii equation (GPE) 



\dtil){yi,t) 



V'(x,t). 
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(In our units h = 1). m is the boson mass, and g quan- 
tifies the interaction strength. We consider a gas in a 
box of size L'^^d ~ 2,3, with periodic boundary condi- 
tions. Diluteness implies that a I, the interparticle 
spacing I — n^^^'^ being determined by the mean density 
n ~ N/L"^. The initial values for the real and imaginary 
parts of the field '0(k, 0) were randomly chosen from a 
Gaussian distribution with width 1/2, centred around 
^n(k,0) cxp{i(f>{k,0)}, where n{k,t) ^ (^/>1"(k, t)V'(k, t)) 
is the occupation number at time t and (/)(k, 0) is a 
random phase angle. Correlation functions, including 
n(k, t), were obtained by averaging over many trajecto- 
ries. We have simulated the time evolution on a grid with 
side length L = NgUs, lattice spacing in d = 2, 3. To 
induce transport from small to large wave numbers, only 
a few modes near k = were chosen to be occupied at the 
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FIG. 1: Phase angle <j){-x.,t) (color scale) at four times during 
a single run of the simulations in d = 2. The white spots 
mark vortex cores where the density falls below 5% of the 



mean density n. Parameters are: g = 2mgal~'^ = 3-10 
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Shown times are: 1. 



t = 26 (top left): Ordered phase shortly after initial prepa- 
ration. 2. t = 820 (top right): After creation of vortex- 
antivortex pairs. 3. t = 6550 (bottom left): During critical 
slowing down of the vortex-antivortex annihilation. 4. t = 10^ 
(bottom right): Low-density vortex-antivortex pairs before fi- 
nal thermalisation. 



initial time with n(k, 0) ^ 1. Such an initial state can 
be prepared, e.g., by Bragg scattering of photons from 
a Bose-Einstein condensate. During the evolution to a 
new equilibrium state, energy is transported to modes 
with higher wave number, and, given appropriate inter- 
action strength, one observes quasistationary momentum 
distributions to develop prior to final thermalisation. 

During the intermediate stage of the equilibration 
quantised vortices appear. To exhibit the vortical 
flow we use the Madelung representation ip{'x,t) = 
y/n{x, t) exp{i(/)(x, <)} of the field in terms of the den- 
sity n(x, t) and a phase angle (/)(x, t). This allows to ex- 
press the particle current j ~ i{^j*'\/^j—ip'Vip*)/2 ~ nw in 
terms of the velocity field v ~ V0. In Fig.[l]wc show the 
phase angle '/'(x, t) during a single run in c? = 2 dimen- 
sions, while in Fig. [2] the evolution in d = 3 is illustrated 
by the near-zero-density points. Many singly quantized 
vortex-antivortex pairs (vortex rings) are formed once the 
initial coherent wave has developed to form shock-wave- 
like fronts. During the later evolution, the pairs mutually 
annihilate at a substantially decreasing rate. After the 
last pair has vanished the gas thermaliscs. In d = 3, de- 
creasingly tangled vortex lines are observed during the 
near-stationary evolution. The initial conditions lead to 
a dense chaotic tangle, as studied also in Refs. pll. [23|. 

Power-law behaviour is found in the momentum dis- 
tributions of particles and energy as we discuss in the 
following. Fig. [3] shows the distribution of occupation 
numbers n{k) = J dild'T-(k) over the radial momentum 
k — |k|, for the four snapshots of the two-dimensional 
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FIG. 2: Snapshots of a single run of the evolution in d = 3 
dimensions. The points show where the density falls below 
5% of the mean density n. Parameters are: g = 8 • 10~* , 
N = 10®, = 128. 1. I = (top left) 2. I = 103 (top right) 
3. I = 410 (bottom left) 4. I = 1640 (bottom right) 

system depicted in Fig. [H on a double-logarithmic scale, 
(dilrf denotes the angular differential.) During the ini- 
tial evolution the mode occupations gradually spread to 
larger wave numbers. As soon as vortex-antivortex pairs 
appear, a power-law regime n{k) = is observed. Dur- 
ing the final stage of the vortex-bearing phase two dis- 
tinct power laws develop which are found to be in excel- 
lent agreement with the analytical prediction of Ref. Q . 
While in the ultraviolet the exponent = d = 2 ex- 
hibits weak wave turbulence, in the infrared the exponent 
confirms the result C = c? + 2 = 4, see Eg. ([2]), corrobo- 
rating results for a relativistic model in [j] . Note that in 
d = 2, the weak-turbulence exponent (^^^ = 2 is identi- 
cal to that in thermal equilibrium in the Rayleigh- Jeans 
regime, n{k) ^ T/uj ^ T/k"^ [2^. The momentum dis- 
tributions during this final phase are shown again, for 
d = 2 and d = 3, in Figs. [4] and [5] (filled black cir- 
cles), respectively, confirming the analytical prediction 
of Eq. ([2|) also for d = 'i. The scaling n{k) ~ fc^^ in 
the ultraviolet reflects that the corresponding modes are 
already thermalised. We remark that, for weaker inter- 
action strengths, the weak-wave-turbulence scaling with 
C,^^ = d = 3 is also seen, at intermediate times. 

To understand these findings in the context of QT 
we analyse kinetic-energy spectra as in [isj . The to- 
tal kinetic energy E'kin = / d'^a; (|V'(/'(x, t)p)/2 can be 
split, iJkin = Ey + Eq, into a 'classical' part = 
J d'^x {\y/n-v\'^) /2 and a 'quantum-pressure' component 
Eq = J d'^a; (|V-y/np)/2. The radial energy spectra for 
these fractions involve the Fourier transform of the gen- 
eralised velocities Wv = y/nrv and Wq = V^/n, 

Es{k) = ^J k'~'dn,{\ws{k)\^), S = v,q. (4) 

Note that the superfluid velocity v — V0 of a single 
vortex diverges as 1 /r with the distance r from the vortex 
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FIG. 3: Single-particle mode occupation numbers as func- 
tions of the radial momentum k = [X^iLi 4sin^(fci/2)]^''^, k = 
27rn/iV,, n = (ni,...,nd), m = -Ns/2, Ns/2, for the four 
different times of the run in d = 2 dimensions shown in Fig. [T] 
Note the double- logarithmic scale. An early development of 
a scaling n{k) ~ is followed by a bimodal scaling with 
n{k) ~ at larger wave numbers. 
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FIG. 4: Single-particle occupation numbers of different frac- 
tions of the system, at time t — 10^ of the d = 2 run shown in 
Fig. [3l The black points are the same as those shown in the 
lower right panel of Fig. [S] Colors distinguish the fractions 
rii (red circles), i.e., the divergence-free part of the flow field 
Wv — y/nv, ric (filled blue squares), i.e., the solenoidal part 
of Wv, and the quantum-pressure partriq (open grey squares). 
A scaling with fc"'* ''^ corresponds to a 5/3 law for the kinetic 
energy in d = 2. = 2n/^ marks the scale corresponding to 
the healing length ^ = l/^/2mgn 



core. Hence, the Fourier transform of v is ill-defined 
and requires regularisation of |v(x)| in the vicinity of the 
core. In Wv, this is naturally achieved by the factor yn 
which rises linearly in r for small r. Following Ref. [l3| 
the regularised velocity Wv is furthermore decomposed 
into 'incompressible' (divergence free) and 'compressible' 
(solenoidal) parts, Wv = Wi -I- Wc, with V • W; = 0, V x 
Wc = 0, to distinguish vortical superfluid and rotationless 
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FIG. 5: Same as in Fig. |4l for the run in d = 3 dimensions 
shown in Fig. [2] at the time t — 1640 of the situation in its 
lower right panel. A scaling with k~^'^^ corresponds to a 5/3 
law for the kinetic energy in d = 3. 

motion of the fluid. 

For comparison of the kinetic-energy spectrum with 
the single-particle spectra n(k), we determine occupation 
numbers corresponding to the different energy fractions 
as ns{k) = k~'^~^Es{k), 6 € {i, c, q}. In the cases we con- 
sider, the resulting spectra ni(fc), nc(fc), and nq(fc) add 
up to the single-particle spectrum n{k) discussed before. 
The components and their sum are shown seperately, for 
d — 2 and d = 3, in Figs. 2] and [51 respectively. Red cir- 
cles denote rii, filled blue squares show the dependence 
of ric, and open grey squares that of riq. Qualitatively, 
the results are the same for d ~ 2 and d ~ 3. Excita- 
tions with large wave numbers are thermally distributed. 
In this regime, the spectrum n{k) k~^ is dominated 
by the compressible and quantum-pressure components. 
For smaller wave numbers the scaling changes to n{k) ^ 
k~'^~'^, being dominated by the velocity Wv. Moreover, 
we find that it is this decomposition into rij and ric which, 
for intermediate wave numbers, allows the incompressible 
part of the energy to develop a Kolmogorov-like scaling 
~ ]^-5/3-d-i a^Qye the scale ki ~ 2tt/1 determined by the 
mean distance I between vortex cores. While the scaling 
of the sum of these components is predicted by the field- 
theoretic analysis in Ref. a rising compressible part 
allows the incompressible contribution to deviate from 
the IR power law and to develop the observed scaling. 
Towards the IR limit, the compressible part becomes too 



weak such that the scaling of ni goes over to C = + 2. 

One can show that the analytically predicted infrared 
power law n{k) ^ fc~^ in d = 2 is consistent with a finite 
density of independent vortices and antivortices. The IR 
scaling ~ k~'^~^ of the compressible (blue) component 
sug gest s an interpretation in terms of acoustic turbulence 
[20I |22|. [25} . Our results show that this component sur- 
vives for a limited period beyond the time when all vorti- 
cal excitations have mutually annihilated. The observed 
scaling corroborates the numerical findings of Ref. [27[ . 
whereby we refrain from sharing the interpretation of the 
power law. 

In summary, our results show a distinct power-law be- 
haviour k^'' of the single-particle momentum spectrum 
n{k) as well as of different components of the kinetic- 
energy distribution over the radial wave number k. Scal- 
ing exponents C of n(k) corroborate the analytical predic- 
tions of Ref. p] . Our findings suggest that local field ex- 
pectation values and short- to intermediate-range coher- 
ence, including topological excitations, are at the basis of 
the infrared power laws predicted within nonperturbative 
dynamical field theory Q . For the chosen generic 
initial conditions, excitations on the top of this classical 
field background support an interpretation in terms of 
acoustic wave turbulence. The connection of these phe- 
nomena with ab-initio dynamical field theoretic methods 
points a way to unified analytical studies of turbulence. 
Moreover, it provides hints of how the proposed dynami- 
cal critical phenomena in relativistic systems [1, [^, Q can 
be realized in nature. 

Experimental studies of universal phenomena in 
nonequilibrium dynamics of ultracold atoms have great 
potential since universal effects do not depend signifi- 
cantly on initial conditions and details of the system. 
The study of turbulence in ultracold gases may have great 
impact on many other fields of physics. Prominent ex- 
amples are strongly correlated nuclear matter produced 
in heavy-ion collisions and early-universe cosmology. 
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